This is a potpourri PR, I will try to move to smaller and more clean PRs in the future. The main changes are:

• Rename ReferenceModel class to ProjectionPredictive
• Use SLSQP optimization method instead of Powell. This takes half the time and the marginals of the projected posteriors look "smooth" and "Gaussian-like", previously it look suspiciously weird (often multimodal or with high peaks and very heavy tails). I will keep exploring other minimization options in the future.
• rename loo_compare to plot_compare we can add an argument in the future if we make the comparison in terms of something else than LOO. I also refactored this function, to reduce the internal dispersion of the code.
• small changes into the docs.